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ABSTRACT 

We explore the quantified morphology of atomic hydrogen (Hi) disks in the Virgo 
cluster. These galaxies display a wealth of phenomena in their Hi morphology, e.g., 
tails, truncation and warps. These morphological disturbances are related to the ram- 
pressure stripping and tidal interaction that galaxies undergo in this dense cluster 
environment. To quantify the morphological transformation of the H i disks, we com- 
pute the morphological parameters of Concentration- Asymmetry-Smoothness, Gini 
and M20 and our own Gm for 51 galaxies in 48 Hi column density maps from the 
VLA Imaging of Virgo spirals in Atomic gas (VIVA) project. 

Some morphological phenomena can be identified in this space of relatively low 
resolution H i data. Truncation of the H i disk can be cleanly identified via the Concen- 
tration parameter (C<1) and Concentration can also be used to identify Hi deficient 
disks (1<C<5). Tidal interaction is typically identified using combinations of these 
morphological parameters, applied to (optical) images of galaxies. We find that some 
selection criteria (Gini-Af2o, Asymmetry, and a modified Concentration- A/20) are still 
applicable for the coarse (~15" FWHM) VIVA Hi data. We note that Asymmetry 
is strongly affected by the choice for the center of these galaxies. The phenomena of 
tidal tails can be reasonably well identified using the Gini-M2o criterion (60 % of 
galaxies with tails identified but with as many contaminants). 

Ram-pressure does move Hi disks into and out of most of our interaction cri- 
teria: the ram-pressure sequence identified by Vollmer et al. (2009) tracks into and 
out of some of these criteria (Asymmetry based and the Gini-M2o selections, but not 
the Concentration-M2o or the Gnbased ones). Therefore, future searches for interac- 
tion using Hi morphologies should take ram-pressure into account as a mechanism 
to disturb Hi disks enough to make them appear as gravitationally interacting. One 
mechanism would be to remove all the Hi deficient (C<5) disks from the sample, as 
these have undergone more than one H i removal mechanism. 
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1 INTRODUCTION 

The Virgo cluster represents an ideal laboratory for the 
study of galaxy evolution in a cluster environment. It is 
one of the nearest galaxy clusters (D~17 Mpc), a rela- 
tively populous system (Abell richness class I), consisting 
of some two thousand catalog ued members i n the Virgo 
Cluster Catalogue (VCC, Bing geU et al] 1 19851 ). The Virgo 
cluster is spiral-rich and dynamically young, i.e., made of 



galaxy subgroups which are falling into the main cluster. 
The most prominent groups are those associated with M86 
to the west and M49 to the south. A particularity of the 
Virgo cluster is that the distribution of its hot intracluster 
medium is strongly peaked on the central cD galaxy M87 
(Bohringor et al. 1994), and overall highly sub-structured. 
The Virgo cluster is therefore ideal to study the effects of 
both gravitational interaction and interaction with the in- 
tracluster medium on the atomic gas (Hi) disks of spiral 
members. 
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Morphological signs of tidal interaction are usually de- 
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tected in the restframe ultraviolet or blue-filter optical im- 
ages of the stellar disk. A series of morphological param- 
eters, originally developed for galaxy classification, have 
been used extensively in recent years to identify disturbed 
galaxies, and to infer the merger fraction s and rates over 
cosmological time sca les fcotz et al J ^200 ^. l2008l: IConselic^ 



20031: IConselice et al.l 12 003. 2008, 2009; 



Yan et al 



Bundv et al.ll2005l: 



2005; Ravindra nath et all |2006| ; iTruiillo et all 



20071 ; I Jogee et al. 2009 1 . The observational benefits are that 



the Hubble Space Telescope observations of high-redshift 
galaxies and ground-based observations of nearby galaxies 
(e.g., SDSS) probe the same wavelength, physical structure 
and volumes. However, in the nearby Universe, the 21 cm 
perspective on galaxies is set to catch up to the optical sur- 
veys in terms of numbers of galaxies studied. Currently un- 
der construction are two new radio observatories: the South 
African Karoo Array Telescope (Me erKAT; Jonas 2 007. ; 
iBooth et al]|2009l : Ide Blok et"al]|2009(), and the Australian 
SKA Pathfi nd er (A SKAP; Ijohnstoni I2OO7I : IJohnston erahl 
I2007l . l2008al lbl. l2009l ). Together with refurbished ob servato- 
ries s uch as the Extended Very Large Array (EVLA: iNapieJ 
I2OO6I ) and the APERture T il e In Focus instrument (APER- 
TIF; IVerheiien et~all I2OO8I : lOosterloo erahl l2009h on the 
Westerbork Synthesis Radio Telescope (WSRT), detailed H I 
maps of thousands of galaxies will become available. The 
unprecedented volume of data implies that the data reduc- 
tion and subsequent morphological classification can only be 

done through automated procedures. 

I n the pr evious papers in this series l|Holwerda et al.l 

I2OO9I . l2011bl ldlfd a), we have shown that the morphological 
classifiers that are in use to identify mergers in the optical 
can be extremely useful for the H I perspective. The H I ob- 
servations are not as high-resolution as the optical disk but 
this is compensated for by the larger extent of the atomic 
gas disk and the relative sensitivity of gas to an interaction - 
for which much anecdotal evide nce already existed i n "The 
H I Rogues Galaxies" catalogue (|Hibbard et al.ll200ll fl. The 
sensitivity of H I surveys to gas-rich and minor mergers will 
provide the calibration for higher redshifts where these types 
of mergers are expected to be the dominant type of galaxy 
inter action jLotz et al.l 2010a |3) ■ The Square Kilometre Ar- 
ray l|Carilli fc Rawlingj 12004 ) will subsequently be able to 
resolve Hi disks at these higher redshifts. 

The Hi disks of spiral galaxies are also extremely sen- 
sitive to ram pressure by the intergalactic or intracluster 
medium. The Virgo cluster is an ideal laboratory to test 
the effects of both gravitational interaction and intraclus- 
ter medium stripping on the quantified morphology of H I 
disks. There are detailed models of the gas stripping se- 
quence (|Vollmeij|2009tV and a uniform H I data set from the 
VLA Imaging of Vi rgo spirals in Atomic gas (VIVA) project 
l|Chung et al.(l2009l ). with detailed notes on individual galax- 
ies identifying morphological phenomena. 

In this paper, we explore the distribution of the H I mor- 
phological parameters fro m the VIVA column density maps 
in relation to estimates bv lChung et"al] (|2007l ). IChung et al.l 
l|2009[) or IVoUmerl l|2009[) of the star-formation rate. Hi de- 
ficiency, stripping and gravitational harassment. First we 
discuss data and morphological parameters (32]and|3]), but 
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readers already familiar with these can skip to results and 
conclusions (Sj4]and[5]). 



2 VIVA DATA 



The VIVA project (jChung et al.l |2009| ) has compiled the 
largest high-resolution imaging H I database of the Virgo 
cluster to date. They have published 48 VLA H I data cubes 
and maps with 51 Virgo cluster spiral galaxies. These are 
some 40 new observations with the remainder from repro- 
cessed archival data. We obtained the zero order moment 
maps from http : / /www . astro . yale . edu/ viva/ and con- 
verted to H l column density maps using the expressions in 
IWaher et all (|2008a|). The spatial resolution of VIVA data 
is approximately 15" (or ~1.2 kpc at the distance of Virgo), 
and the typical sensitivity is ~ 3 — 5 x 10^^ cm^'^ in 3a 
per 10 km/s channel. VIVA data show that (i) the galaxies 
located in the core of the cluster predominantly have trun- 
cated H I disks with some gas displaced from the galactic disk 
and (ii) some 7 out of t he 50 sample g alaxies (~10%) show 
long, one-sided Hi tails l|Chung et al.l l2007l. These tails oc- 
cur in galaxies located at intermediate projected distances 
between 0.5 and 1 Mpc fr om the cluste r cent re. From the 
sample listed in Table 1 of IChung et ahl l|2009l ). we did not 
use IC 3418 as it is a non-detection in the VIVA survey and 
Holmberg VII because it does not have enough continuous 
Hi fiux to meaningfully compute morphological parameters. 



3 MORPHOLOGICAL PARAMETERS 

Two morphological parameter schemes are now in 
wide use: Co ncentration-Asymm etry-Smoothness (CAS, 
IConselicd I2OO3I ) and Gini-M2o l|Lotz et al.l 120041 ) . In an 
image with n pixels with intensities at pixel position 

CAS the morphological parameters are defined as: 



C = 5 



A = 



S».,|J(i,j)-Ji8o(i,j)| 
S,..|/(i,j)l 

^^,Mhj)-Is{hj)\ 

s,,.|/(i,j)| 



(1) 

(2) 
(3) 



where r% is the radius which includes that percentage of the 
intensity of the object, Iiso{i,j) is the value of the pixel in 
the rotated image and Is(i,j) is the same pixel in the image 
after smoothing. 

Gini and M20 are defined as: 



G 



M20 = log 



In{n - 1) 



E.(2i-n-l)|J,l (4) 
for which E*^7i < 0.2 Itot is true. (5) 



with Mtot = SMi = T,Ii[{xi-Xc)^ + {yi~yc)^], where {xc,yc) 
is the galaxy's central position, and U is the intensity of pixel 
i in an flux-ordered list. Pixel k is the pixel marking the 20% 
point in this list. / is the mean pixel value. The Gini param- 
eter is an indicator of inequality in the distribution of pixel 
values; Gini=l is perfect inequality with one pixel contain- 
ing all the flux from an object and Gini=0 is perfect equality 
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with each pixel containing an equal fraction of the flux. Val- 
ues in between indicate how concentrated the flux from this 
object is in a few bright areas. The M20 parameter is an 
indication whether or not the brightest parts of the object 
are close to the center or spread throughout the object. 

In addition to these well-established parameters, we 
compute the Gini coefficient of the second or der moment Ij x 
\(xi — Xc)^ + (yi — yc)^], which we called Gm l|Holwerda et al.l 
l201ld ). For all these parameters, we compute the uncer- 
tainty by varying the input centre the resolution and ran- 
domizing the pixel values. The exception is the Gini param- 
eter as this does not rely on the centre as input. For the 
Gini error, we jacknife the pixel selection. For a more de- 
tailed description of these paramet ers, we refer the reader to 
the previous papers in this series (iHolwerda et al.ll2011bllH ) 
and the origina l papers introducing these jConselice 200c ; 
iLotz et al]|2004l '). We compute these parameters for the area 
defined by the contour l evel of 1 x 10^^ cm~^ (the la level). 

In iHolwerda et al.l l|2011cl ). we identified several parts 
of parameter space where tidally interacting galaxies reside 
for the WHISP sample. To identify the interacting galaxies, 
we used a subsample of 154 spirals and dwarf galaxies in 
WHISP, which have independent visu al estimates of interac- 
tion fr om either [Swaters et al.l (|2002l ) or lNoordermeer et al.l 
(|2005l ). The criteria from the literature and our previous pa- 
per are identified in Figures[2]and[3]with dashed and dotted 
lines respectively. Criteria from the literature are: 

A > 0.38 (6) 

ft-om lConselic3 l|2003l ) , marked by the dashed lines in Figures 
[land El panels IV, V and VI. 

G > -0.115 X M20 0.384 (7) 

from lLotz et~al] l|2004l ). marked by the dashed lines in panel 
II of Figures [3 and [31 

G > -0.4 X A 0.66 or A > 0.4 (8) 

iLotz et akl l|2010al ). marked by the dashed lines in panel IV 
of Figures [Hand [3 

We defined three criteria based on the WHISP H I mor- 
phologies: 



Gm > 0.6, 



(9) 



makred by the vertical dotted line in Figures [2l and [3l panels 
I, HI, VI and X, 

A > -0.2 X M20 + 0.25, (10) 
marked by the dotted line in Figures [2l and [3l panel V, and. 



G > -5 X M20 -f 3, 



(11) 



which is the dotted lines in Figures [21 and [3 panel IX. 

These criteria select a merger from its H I morp hol- 
ogy for different time scales. In Holwcrda et al.l ([201l3)) we 
quantified the typical time scale for a merger to be selected 
by each of these criteria using Hi maps of me rger simu- 
lations (originally presented in ICox et al.l [2006 a bV These 
can be compared to merger simulations of an isolated disk, 
evolving passively. The comparison shows how long a merger 
is selected as well as the level of contamination by passive 
disks. At the resolution of the WHISP sample, interactions 
were identifiable as such by their H l morphology for about 
a Gyr. This is a similar time scale as the one for optical 



morphological identifications and better than the close pair 
method. Our goal is to explore if the above criteria could 
be applied to Hi data in a denser environment, where many 
other phenomena play a role in the H I morphology as well 
as to explore how many of those other phenomena leave a 
signature in the H I morphology. 



4 RESULTS: Hi MORPHOLOGY 

Figures [2l and [3l show the spread of morphological parame- 
ters for the entire V IVA sample with different morphological 
classifications from IChung et all (|2007l . |2009| '1. Table [1] lists 
all the values. Due to the combination of distance and res- 
olution, the VIVA sample has the coarsest physical scales 
we have applied these m orphological parameters to d ate. 
In IHolwerda et al.l (|2011bl ') and Holwerda et al.l (|2011dl ') we 
showed that these parameters are relatively invariant with 
distances out to D~20 Mpc and a resolution of 6", so the 
VIVA sample is effectively at the limit of this technique with 
current observatories. 

We us ed the optical posit ions (J2000) of the galaxies 
reported in IChung et al.) (|2009l ') as the centre of the galaxy 
{xc,yc) to compute the parameters that depend on a central 
position (all except the Gini paramter). Our motivations are 
that the stellar disk marks the centre of the gravitational 
potential of a galax y and kinematic and stellar centres gen- 
erally align well fsee lWalter et al]l2008bl : lTrachternach et al.l 
(iooi). Optical positions of the galaxies are a l ikely starting 
point for many of the future H I surveys and IChung et al.l 
(|2009l ) reported many tidal Hi tails and displacements, 
which could be quantified using our morphological parame- 
ters. Thus, it is entirely possible that some of morphological 
values (e.g.. Asymmetry) are extreme since the centre is not 
the barycentre of the fiux, as is common for optical morpho- 
logi cal measurem e nts. 

IChung et al.1 tOO^ supply extensive notes on the H I 
morphology and the nature of the H I stripping of individual 
galaxies. We have summarized these evaluations in Table [5l 
I Chung et al.. (2009 ) note which galaxies are currently under- 
going or have recently undergone stripping, as well as the 
likely origin of this stripping; ram pressure by the intraclus- 
ter medium or tidal interaction with a close companion. In 
addition, they report a parameter for Hi deficiency [defni', 
the galaxy's Hi surface density compared to the value typ- 
ical for spirals), a measure of Hi truncation (Dhi/Db; the 
ratio of the Hi diameter over the optical, B-band one), and 
additional notes for indicators of warps in the H I morphol- 
ogy and velocity maps. Finally, IChung et al.l (|2009l ') also list 
the star-formatio n properties of their sample, a s catalogued 
and classified bv lKoopmann fc KennevI ()2004al lbll: the level 
of activity and whether it was recently cut-off (truncated 
star- formation). We have plot the spread of our morpho- 
logical parameters with each of these estimates marked to 
explore if there are specific parts of parameter space that 
single out an Hi disk characteristic, similar to the separa- 
tion of interaction and non-interacting galaxies we found in 
the WHISP database. Figures [T] [2l and[3l show the spread 
in parameters each with certain phenomena marked, so that 
they can be identified in the morphological parameter space. 
We discuss several Hi morphological phenomena noted by 
Chung et al. below. 
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Figure 1. The relation between Gini and Concentration for Hi truncated galaxies {Dhi ^ Db, left panel) and Hi deficient galaxies 
(right panel). Most of the VIVA sample disks are Hi deficient (defui < 1) and we mark t hose that are extremely so (defui < 0.5); 
galaxies with half a dex or more below the typical H I surfac e density for their type of spiral l lHavnes fc Giovanelli|[l984h . Typical values 
for concentration are 5 and above for Hi disks in WHISP jH olwerda et al.ll201l3 lah. Therefore the majority of the H I disks in VIVA 
are more than typically centrally concentrated; a low concentration value (C< 5) is a good sign of Hi deficiency and, in addition, an 
extremely low concentration value (C<1) identifies half of the truncated Hi disks. 



4.1 Truncation and Deficiency 

Figure [1] shows Gini versus Concentration for H I deficient 
and truncated Hi dislts. Botli parameters are sensitive to 
iiow evenly tlie flux is distributed over the area of tlie disk 
and tlierefore could be sensitive to botli tliese phenomena. 
To identify Hi truncation, one could use tlie Concentration 
of the Hi emission. We have found tliat values around or 
below C=l seem to be very indicative of truncation; the 
radii containing 20 and 80 percent of the Hi flux are very 
close together (Fig. [T] left panel). Given that typical val- 
ues for concentration i n Hi are much higher (typically C^5, 
iHolwerda et al.l|201ld f3). this may be a good tracer of trun- 
cation in large Hi surveys. A value of Concentration be- 
low unity selects half of the truncated galaxies in the VIVA 
sample very cleanly (12 out of 24 truncated disks with no 
contamination from non-truncated disks). 

The H I deficiency is defined by the surface density com- 
par ed to the average su rface density for spirals (see equation 
5 in lChung et al.ir2009l ). The majority of VIVA galaxies are 
Hi deficient (44 galaxies with defni < 1). The low values 
of concentration (C < 5) appear to be indicative of this H I 
deficiency. In VIVA, those galaxies with extreme Hi defi- 
ciency (defHi < 0.5) are not much more concentrated than 
moderately deficient ones (Figure [1] right panel), at least 
at this resolution. We note that a C^5 criterion will select 
some 80 % of the H I deficient galaxies in the VIVA sample 
but their number is contaminated (17% of this selection is 
not Hi deficient). In the upcoming large surveys this makes 
Concentration a quick assessor on whether or not a spiral 
has a typical Hi content or not. 



4.2 Interaction; H I Tails, and Ram-pressure 

Nearly all the spirals in the Virgo cluster are undergoing, 
or have recently undergone, some type of interaction. Dis- 
cerning between the dominant mechanism, t idal or ram- 
pressu re, was a major fo cus of the work by IChung et al.l 
(|2009[ ) and lVollmeJ (|2009l ). Several interaction related mor- 
ph ological phen o mena (e.g., warps and Hi tails) are noted 
bv lChung et"al] l|2007l ). 

Table [S] lists many galaxies with H I tails in the VIVA 
data and several with H I warps. One would expect that high 
values for the Asymmetry would identify Hi tails (Figure 
21). However, because we used the optical centres listed by 
Chung et al.l (|2009, ) , Asymmetry is already at the maximum 
value for most of the VIVA galaxies (Figure [2l panels IV- 
VII). We note however that the Asymmetry may still be a 
good indicator of (tidal) H I tails, but in less extreme or pos- 
sibly better resolved cases. The M2o-Gini selection criterion 
(see eq. O does reasonably separate those galaxies with H I 
tai ls from those which do not ; 60 % of the Hi tails identified 
by IChung et al.1 (|2007l . l2009h are selected by this criterion 
(Figure [21 panel II). However, these are about half the ob- 
jects selected do not have an H I tail. Therefore, the Gini-M2o 
criterion can be used a likelyhood for an H I tail but it will 
suffer from contamination. Warps were identified predomi - 
nantly by their kinematic signature in IChung et al.l (|2009l ) 
and no morphological parameter or combination seems to 
be able to identify them in the VIVA data. 

Figure O shows the same spread of morphological pa- 
rameters as Figure [2l with the estimate we gleaned from 
IChung et al.l l|2009l ) whether or not the galaxy is stripped. 
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Figure 2. The spread of morphological parameters. Galaxies with an Hi tail, as identified in lChung et al ] l|2007l ) and lChung et al ] l|2009l ') 
are marked with solid black points. The (optical) merger selection criteria from the literature are marked with dashed line s in panel II 
(equation O , panel IV (equations [6] and (Sjl , and V and VI (equation (6)1 . Our selection criteria from iHolwerda et al.l l|201ld ) are marked 
with dotted lines; the Gm criterion in panels I, III, VI and X (equation O , the j4-M2o criterion in panel V mOII and the C-M20 criterion 
in panel IX (equation [TT} . 



undergoing ram-pressure, or undergoing tidal interaction. 
The seq uence of stripping in the Virgo cluster according to 
IVollmeii (j2009) is indicated by the solid lines. This sequence 
of six Virgo galaxies does trace a track into and out of cer- 
tain interaction criteria defined previously in the literature 
and the third paper in this series. 

The extreme values for Asymmetry drive almost all of 
the VIVA galaxies into the Asymmetry-based merger selec- 



tion space; the Asymmetry criterior i from | Conselicel (l2003h . 
the Asymmetry-Gini criterion from'Lotz et al.' (^OlOa") and 
our Asymmetry-M2o criterion from Holwcrda et al. (2 01l3 ) 
(equations miHl and[lO]respectively). The extreme values for 
Concentration -linked to H I deficiency and truncation- sim- 
ilarly m ake the Concentration- M2o selection criterion (equa- 
tion [m iHolwerda et~al]l2011c[ l not applicable for merger se- 
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Figure 3. The spread of morphological parameters for gala xies underKo ing tidal interaction (black symbols), Hi stripping (open circles) 
and ongoing ram-pressure (marked with an x), according to lChung et al. (2009), see TablefS] The (optical) merger selection criteria from 
the literature are ma rked with dashed lines in panel II (equation (T)! , panel IV (equations |6] and [8]l , and V and VI (equation |6]l . Our 
selection criteria from lHoIwerda et al. I ll201ld'l are marked with dotted lines; the Gm criterion in panels I, III, VI and X (equation |9ll , the 
A-M20 criterio n in pan e l V II10|I and the C-M20 criterion in panel IX (equation II 111 . The solid lines connect the ram-pressure sequence 
determined by IVoUmed l l200g|'l for a select few Virgo cluster galaxies. 



lection. Only the Gini-M2o criterion (eq. [7| appears to sep- 
arate objects with Hi tails moderately as discussed above. 

The Gm criterion (eq. [9| is not contaminated by galax- 
ies undergoing ram-pressure stripping but does not select 
tidal interactions well anymore either, something that we 
attribute to the native VIVA res olution (no visibility times 
found in iHolwerda et al.ll2011dl . for the VIVA resolution). 



Because stripping. Hi deficiency and truncation are inter- 
linked, there is the possibility that the Concentration selec- 
tion criteria also will select those galaxies that are or have 
been stripped. 

Two effects move Hi morphologies into the selection 
criteria: poorer resolution and additional effects such as the 
ram-pressure stripping. Given how the VoUmer sequence 
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tracks into (and out again) of the merger selection crite- 
ria, we can conclude that galaxies undergoing ram-pressure 
stripping will sometimes be selected as interacting in a mor- 
phological selection in future Hi surveys, at least those with 
spatial resolutions similar to VIVA. The addition of some 
galajcies that are undergoing stripping would, in part, ex- 
plain the high merger fraction and ra te we found in the 
WHISP survey l|Holwerda et al.ll2011al '). In iHolwerda et al] 
(|2011dl ). we noted that the lower resolution data of the VIVA 
sample would have some problems selecting interacting sys- 
tems from isolated ones, purely on the basis of their H I 
morphology (Table 8 in the appendix). However, the success 
with which tidal H I tails could be separated is encouraging. 

Therefore, in future medium-resolution H I surveys (e.g., 
WALLABY, Koribalski et al, in preparation), the selection 
of Hi disks with disturbed morphology can be used as a 
first pass in a search for stripping or interacting galaxies 
and a merger fraction and rate will have to be corrected 
for the contamination by stripping galaxies. For example, 
those galaxies that are Hi deficient can be removed before 
a merger rate is determined but that would bias against 
denser cluster environments (which do include more other 
interactions like ram-pressure stripping). A higher resolu- 
tion survey of the Virgo cluster galaxies would resolve the 
population of tidally interacting galaxies. The separation of 
ram-pressure stripped galaxies from tidal interacting ones 
may then be possible without the removal of Hi deficient 



4.3 Star Formation 

We explored the link between the morphological parameter 
space of H I and the star-format ion and truncation of star- 
formation as listed in Table 3 in lChung et al.l (j2009'). Rela- 
tions are weak as the Hi gas is undergoing several other ef- 
fects other than feedback from star-formation but it is inter- 
esting to note that those galaxies with little star-formation 
also have lower values of Gini for Hi (i.e., they are more 
homogeneous) . 



and M20 interaction criterio n (eq. [7] & Figure [2] panel II) 
from lLotz et~all (|2004 l2008h . 

Galaxies that are interacting are still identifie d by some 
of th e morphological c riteria from the literature (jConselicel 
l2003l : iLotz et aL 12004]') and those id entified by us specifi- 
cally for H I in lHolwerda et al] (|2011cl ) (eq.[6]-[ni- However, 
galaxies affected by ram-pressure stripping are ofttimes also 
selected by these criteria. For exa mple, we note t hat the 
ram pressure sequence identified by IVoUmerl (|2009l ) tracks 
into the interaction part of the parameter space for many 
of the morphological selection criteria (solid tracks in Fig- 
ure |3} . Discerning between the origin of the H I stripping 
-ram pressure or tidal- is, however, not possible with these 
morphological parameters, at least at VIVA's spatial reso- 
lution. Thus any selection of interacting systems in a survey 
based purely on the H I morphology parameters will include 
those that are undergoing ram-pressure stripping in a denser 
environment. A straightforward correction would be to ex- 
clude all H I deficient (C < 5) galaxies from a morphological 
selection of interacting galaxies. 
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5 CONCLUDING REMARKS 

In this paper, we present the morphological parameters com- 
monly used in the optical, applied to the Hi column den- 
sity maps of Virgo Cluster galaxies obtained by the VIVA 
program. We find that the resolution of these observations 
(^15" or 1 kpc) is on the coarse side for this type of analy- 
sis but some Hi morphological phenomena are identifiable, 
though with limitations. A critical choice is whether or not 
to use the optical or the Hi centres for these galaxies, one 
that every H I survey will have to make. The choice of centre 
affects measures such as concentration, M20, Gm, and espe- 
cially Asymmetry. Alternatively, one could use the barycen- 
tre of the H I emission itself but the optical centre is most 
likely the centre of the dark matter halo of the galaxy and 
this approach would likely miss Hi displacements like tidal 
tails. We find that of the Hi morphological phenomena. Hi 
deficiency and truncation can be identified easily by a low 
value of the concentration index (Figure[T|. Warps could not 
be identified using the morphology parameters and H I tails 
can be identified moderately well with the combined Gini 



REFERENCES 

Binggeh B., Sandage A., Tammann G. A., 1985, AJ, 90, 
1681 

Bohringer H., Briel U. G., Schwarz R. A., Voges W., Hart- 
ner G., Triimper J., 1994, Nature, 368, 828 

Booth R. S., de Blok W. J. G., Jonas J. L., Fanaroff B., 
2009, ArXiv e-prints/0910.2935 

Bundy K., Ellis R. S., Conselice C. J., 2005, ApJ, 625, 621 

Carilli C. L., Rawlings S., 2004, New Astronomy Review, 
48, 979 

Chung A., van Gorkom J. H., Kenney J. D. P., Growl H., 

VoUmer B., 2009, ArXiv e-prints 
Chung A., van Gorkom J. H., Kenney J. D. P., VoUmer B., 

2007, ApJ, 659, L115 
Conselice C. J., 2003, ApJS, 147, 1 

Conselice C. J., Bershady M. A., Dickinson M., Papovich 

C., 2003, AJ, 126, 1183 
Conselice C. J., Bundy K., U V., Eisenhardt P., Lotz J., 

Newman J., 2008, MNRAS, 383, 1366 
Conselice C. J., Yang C, Bluck A. F. L., 2009, MNRAS, 

361 



8 B. W. Holwerda et al. 



Cox T. J., Dutta S. N., Di Matteo T., Hernquist L., Hop- 
kins P. F., Robertson B., Springel V., 2006a, ApJ, 650, 
791 

Cox T. J., Jonsson P., Primack J. R., Somerville R. S., 
2006b, MNRAS, 373, 1013 

de Blok W. J. C, Jonas J., Fanaroff B., Holwerda B. W., 
Bouchard A., Blyth S., van der Heyden K., Prizkal N., 
2009, in Conference Proceedings of the "Panoramic Radio 
Astronomy: Wide-field 1-2 GHz research on galaxy evolu- 
tion" , June 02 - 05, 2009 Groningen, the Netherlands 

Haynes M. P., GiovaneUi R., 1984, AJ, 89, 758 

Hibbard J. E., van Gorkom J. H., Rupen M. P., Schimi- 
novich D., 2001, in Astronomical Society of the Pacific 
Conference Series, Vol. 240, Gas and Gala:xy Evolution, 
Hibbard J. E., Rupen M., van Gorkom J. H., eds., pp. 
657— h 

Holwerda B. W., de Blok W. J. G., Bouchard A., Blyth 

S., van der Hey den K., Prizkal N., 2009, in Conference 
Proceedings of the "Panoramic Radio Astronomy: Wide- 
field 1-2 GHz research on galaxy evolution" , June 02 - 05, 
2009 Groningen, the Netherlands 
Holwerda B. W., Pirzkal N., de Blok W. J. G., Blyth S.- 
L., Bouchard A., van der Heyden K. J., 2011a, MNRAS, 
submitted 

Holwerda B. W., Pirzkal N., de Blok W. J. G., Blyth S.-L., 
Bouchard A., van der Heyden K. J., Elson E. C, 2011b, 
MNRAS, accepted 

—, 2011c, MNRAS, accepted 

— , 2011d, MNRAS, submitted 

Holwerda B. W., Pirzkal N., de Blok W. J. G., van Driel 

W., 2011e, MNRAS, in preparation 
Jogee S., Miller S. H., Penner K., Skclton R. E., Consclice 

C. J., Somerville R. S., BcU E. F., Zheng X. Z., Rix H., 
Robaina A. R., Barazza F. D., Barden M., Borch A., Beck- 
with S. V. W., Caldwell J. A. R., Peng C. Y., Heymans C, 
Mcintosh D. H., Haufiler B., Jalmkc K., Mciscnheimer K., 
Sanchez S. F., Wisotzki L., Wolf C, Papovich C, 2009, 
ApJ, 697, 1971 

Johnston S., 2007, in From Planets to Dark Energy: the 
Modern Radio Universe. October 1-5 2007, The Univer- 
sity of Manchester, UK. Published online at SISSA, Pro- 
ceedings of Science, p. 6 

Johnston S., Bailes M., Bartel N., Baugh C, Bietenholz 
M., Blake C, Braun R., Brown J., Chatterjee S., Dar- 
ling J., Deller A., Dodson R., Edwards P. G., Ekers R., 
EUingsen S., Fcain I., Gaensler B. M., Havcrkorn M., 
Hobbs G., Hopkins A., Jackson C, James C, Joncas G., 
Kaspi v., Kilborn V., Koribalski B., Kothes R., Landecker 
T. L., Lenc E., Lovell J., Macquart J.-P., Manchester R., 
Matthews D., McClurc- Griffiths N. M., Norris R., Pen U.- 
L., Phillips C, Power C, Prothcroe R., Sadler E., Schmidt 
B., Stairs I., Staveley-Smith L., Stil J., Taylor R., Tingay 
S., Tzioumis A., Walker M., Wall J., WoUeben M., 2007, 
Publications of the Astronomical Society of Australia, 24, 
174 

Johnston S., Fcain I. J., Gupta N., 2009, in Astronom- 
ical Society of the Pacific Conference Series, Vol. 407, 
Astronomical Society of the Pacific Conference Series, 

D. J. Saikia, D. A. Green, Y. Gupta, & T. Venturi, ed., 
pp. 446—1- 

Johnston S., Taylor R., Bailes M., Bartel N., Baugh C, 
Bietenholz M., Blake C, Braun R., Brown J., Chatter- 



jee S., Darhng J., DcUer A., Dodson R., Edwards P., Ek- 
ers R., EUingsen S., Fcain I., Gaensler B., Haverkorn M., 
Hobbs G., Hopkins A., Jackson C, James C, Joncas G., 
Kaspi v., Kilborn V., Koribalski B., Kothes R., Landecker 
T., Lenc A., Lovell J., Macquart J.-P., Manchester R., 
Matthews D., McClure-Griflaths N., Norris R., Pen U.-L., 
Phillips C, Power C, Protheroe R., Sadler E., Schmidt B., 
Stairs I., Staveley-Smith L., Stil J., Tingay S., Tzioumis 
A., Walker M., Wall J., WoUeben M., 2008a, Experimental 
Astronomy, 22, 151 

— , 2008b, Experimental Astronomy, 22, 151 

Jonas J., 2007, in From Planets to Dark Energy: the Mod- 
ern Radio Universe. October 1-5 2007, The University of 
Manchester, UK. Published online at SISSA, Proceedings 
of Science, p. 7 

Koopmann R. A., Kenney J. D. P., 2004a, ApJ, 613, 866 

— , 2004b, ApJ, 613, 851 

Lotz J. M., Davis M., Faber S. M., Guhathakurta P., Gwyn 
S., Huang J., Koo D. C, Le Floc'h E., Lin L., Newman 
J., Noeske K., Papovich C, Willmcr C. N. A., CoU A., 
Conselice C. J., Cooper M., Hopkins A. M., Metevier A., 
Primack J., Rieke G., Weiner B. J., 2008, ApJ, 672, 177 

Lotz J. M., Jonsson P., Cox T. J., Primack J. R., 2010a, 
MNRAS, 404, 590 

— , 2010b, MNRAS, 404, 575 

Lotz J. M., Primack J., Madau P., 2004, A J, 128, 163 

Napier P. J., 2006, in Astronomical Society of the Pa- 
cific Conference Series, Vol. 356, Revealing the Molecular 

Universe: One Antenna is Never Enough, D. C. Backer, 
J. M. Moran, & J. L. Turner, ed., pp. 65— |- 

Noordermeer E., van der Hulst J. M., Sancisi R., Swaters 
R. A., van Albada T. S., 2005, A&A. 442, 137 

Oosterloo T., Verheijen M., van Cappellen W., Bakker L., 
Heald G., Ivashina M., 2009, ArXiv e-prints 

Ravindranath S., Giavalisco M., Ferguson H. C, Consclice 
C, Katz N., Weinberg M., Lotz J., Dickinson M., FaU 
S. M., Mobasher B., Papovich C, 2006, ApJ, 652, 963 

Swaters R. A., van Albada T. S., van der Hulst J. M., 

Sancisi R., 2002, A&A, 390, 829 

Trachternach C, de Blok W. J. G., Walter F., Brinks E., 
Kennicutt R. C, 2008, AJ, 136, 2720 

Trujillo I., Conselice C. J., Bundy K., Cooper M. C, Eisen- 
hardt P., Ellis R. S., 2007, MNRAS, 382, 109 

Verheijen M. A. W., Oosterloo T. A., van Cappellen W. A., 

Bakker L., Ivashina M. V., van der HxUst J. M., 2008, 
in American Institute of Physics Conference Scries, Vol. 
1035, The Evolution of Galaxies Through the Neutral Hy- 
drogen Window, R. Minchin & E. Momjian, ed., pp. 265- 
271 

VoUmer B., 2009, A&A, 502, 427 

Walter F., Brinks E., de Blok W. J. G., Bigiel F., Kennicutt 
R. C, Thornley M. D., Leroy A., 2008a, AJ, 136, 2563 

— , 2008b, AJ, 136, 2563 

Yan H., Dickinson M., Stern D., Eiscnhardt P. R. M., 
Chary R.-R., Giavalisco M., Ferguson H. C, Casertano 
S., Conselice C. J., Papovich C, Reach W. T., Grogin N., 
Moustakas L. A., Ouchi M., 2005, ApJ, 634, 109 



Hi Morphology in Virgo 9 



i> V ' v.. 


V 


stup 




tidal 


warp 


tail 




4064 




/ 




/ 








4189 


89 














4192 


92 








/ 






4216 


167 














4222 


187 








/ 






4254 


307 






/ 




/ 




4293 


460 


/ 




/ 








4294 


465 






/ 




/ 


a 


4298 


483 






/ 








4299 


491 


/ 


/ 


/ 




a 




4302 


497 


/ 


/ 






/ 


a 


4321 


596 






/ 








4330 


630 


/ 


/ 






/ 


a, I 


4351 


692 




/ 


/ 








4380 


792 














4383 


801 






/ 


/ 






4388 


836 


/ 


/ 






/ 


II 


4394 


857 














4396 


865 


/ 


/ 






/ 


a 


4405 


874 


/ 












4402 


873 


/ 


/ 






/ 




IC3355 


945 






/ 








4419 


958 


/ 












4424 


979 


/ 




/ 




/ 


a 


4450 


1110 














IC3392 


1126 


/ 












4457 


1145 














4501 


1401 




/ 








III 


4522 


1516 


/ 


/ 








IV 


4532 


1554 






/ 


/ 


/ 




4535 


1555 




/ 










4533 


1557 






/ 




/ 




4536 


1562 






/ 








4548 


1615 














4561 








/ 








4567 


1673 






/ 




/ 




4568 


1676 






/ 


/ 






4569 


1690 


/ 


/ 








V 


4579 


1727 














4580 


1730 


/ 


/ 










4606 


1859 


/ 




/ 








4607 


1868 














4651 








/ 


/ 


/ 




4654 


1987 


/ 


/ 






/ 


a 


4689 


2058 
2062 






/ 








4694 


2066 






/ 








4698 


2070 




/ 


/ 




/ 




4713 










/ 






4772 








/ 


/ 






4808 








/ 
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Table 1. The morphological values and errors of the VIVA Hi column density maps. 
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